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We show that the NP-hard quadratic unconstrained binary optimization (QUBO) problem on a graph 
G can be solved using an adiabatic quantum computer that implements an Ising spin- 1/2 Hamiltonian, by 
reduction through minor-embedding of G in the quantum hardware graph U . There are two components to 
5^ i this reduction: embedding and parameter setting. The embedding problem is to find a minor-embedding Gemb 

■ of a graph G in U, which is a subgraph of U such that G can be obtained from Gemb by contracting edges. The 

parameter setting problem is to determine the corresponding parameters, qubit biases and coupler strengths, 
of the embedded Ising Hamiltonian. In this paper, we focus on the parameter setting problem. As an example, 
we demonstrate the embedded Ising Hamiltonian for solving the maximum independent set (MIS) problem via 
adiabatic quantum computation (AQC) using an Ising spin- 1/2 system. We close by discussing several related 
algorithmic problems that need to be investigated in order to facilitate the design of adiabatic algorithms and 
QQ I AQC architectures. 

00 ' 

, 1 Introduction 

O . 

OO , Adiabatic Quantum Computation (AQC) was proposed by Farhi et al. |[T5l[T4l in 2000. The AQC model is based 
on the adiabatic theorem (see, e.g. ||22]| ). Informally, the theorem says that if we take a quantum system whose 
Hamiltonian "slowly" changes from T^init (initial Hamiltonian) to T^finai (final Hamiltonian), then if we start 
^ ■ with the system in the groundstate (eigenvector corresponding to the lowest eigenvalue) of T-i\„\t, then at the 
end of the evolution the system will be "predominantly" in the ground state of Tifinai- The theorem is used to 
construct adiabatic algorithms for optimization problems in the following way: The initial Hamiltonian T^init 
is designed such that the system can be readily initialized into its known groundstate, while the groundstate of 
the final Hamiltonian T^finai encodes the answer to the desired optimization problem. The complete (or system) 
Hamiltonian at a time t is then given by 

n{t) = fl-s(4)Vinit + s(;^)7ifinal 



O 



for t € [0, T] where s increases monotonically from s(0) = to s(l) = 1 and T is the total running time of the 
algorithm. If T is large enough, which is determined by the minimum spectral gap (the difference between the 
two lowest energy levels) of the system Hamiltonian, the adiabatic theorem guarantees the state at time t will be 
the groundstate of 7Y(t), leading to the solution, the ground state of 7Y(r) = Hfinai- 
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It is believed that AQC is advantageous over standard (the gate model) quantum computation in that it is 
more robust against environmental noise fTTl [3l l4ll. In 2004, D-Wave Systems Inc. undertook the endeavor 
to build an adiabatic quantum computer for solving NP-hard problems. Kaminsky and Lloyd [17] proposed a 
scalable architecture for AQC of NP-hard problems. Note that AQC is polynomially equivalent to the standard 
quantum computation, and therefore AQC is universal. Farhi et al. in their original paper ifTSi showed that 
AQC can be efficiently simulated by standard quantum computers. In 2004, Aharonov et al fT] proved the more 
difficult converse statement that standard quantum computation can be efficiently simulated by AQC. In this 
paper, however, we will focus on a subclass of Hamiltonians, known as Ising models in a transverse field, that 
are NP-hard but not universal for quantum computation. This subclass of Hamiltonians has been implemented 
by D-Wave Systems Inc. 

D-Wave's quantum hardware architecture can be viewed as an undirected graph U with weighted vertices 
and weighted edges. See Figure[I]for an example. Denote the vertex set of U by V(C/) and the edge set of U by 




Figure 1: An example hardware graph: a 4 
nearest neighbors. 



4 extended grid. Each qubit is coupled with its nearest and next- 



E({7). Each vertex i G V({7) corresponds to a qubit, and each edge ij € E(C/) corresponds to a coupler between 
qubit i and qubit j. In the following, we will use qubit and vertex, and coupler and edge interchangeably when 
there is no confusion. There are two weights, hi (called the qubit bias) and Aj (called the tunneling amplitude), 
associated with each qubit i. There is a weight Jij (called the coupler strength) associated with each coupler ij. 
In general, these weights are functions of time, i.e., they vary over the time, e.g., hi{t). 
The Hamiltonian of the Ising model in a transverse field thus implemented is: 



Hit) 



h,{t)a!+ E 

j6V(G) ij€E{G) i6V(G) 

with af = / / (g) . . . cj^ . . . (X) / (the is in the ith position), similarly for af and crfaj, where I is the 
2x2 identity matrix, while and are the Pauli matrices given by 



1 
-1 



and cr^' 



1 

1 



In general, the transverse field (i.e. ^jgv(G) ^i^f^ encodes the initial Hamiltonian 7^init> while the other two 
terms encode the final Hamiltonian: 



Wfinal = E + E '^ij'^i'^j- 



(1) 
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One can show (see Appendix) that the eigenvalues and the corresponding eigenstates of Wfinai are encoded in the 
following energy function of Ising Model: 



£'(si , . . . , Sn 



) — ^ ^ hi Si -\- ^ ^ JijSiSj 

ieV(G) ueE(G) 



(2) 



where Sj € {—1, +1}, called a spin, and hi, Jij € M. In particular, the smallest eigenvalue of Hf\na\ correspond- 
ing to the minimum of £, and arg minS corresponds to its eigenvector (called ground state) of T^finai- When 
there is no confusion, we use the energy function of Ising model and Ising Hamiltonian interchangeably. Here- 
after, we refer the problem of finding the minimum energy of the Ising Model or equivalently the groundstate of 
Ising Hamiltonian as the Ising problem. It can be shown that the Ising problem is equivalent to the problem of 
Quadratic Unconstrained Boolean Optimization (QUBO)(see 111 and references therein), which has been shown 
to be a common model for a wide variety of discrete optimization problems. More specifically, finding the mini- 
mum of (J in (l2ll is equivalent to finding the maximum of the following function (which is also known as quadratic 
pseudo-Boolean function LSJ) of QUBO on the same graph: 



where Xi G {0, 1}, Cj, Jij G M. The correspondence between the parameters, /is and cs will be shown in Section|2l 
Therefore, given an Ising/QUBO problem on graph G, one can thus solve the problem on an adiabatic quantum 
computer (using an Ising spin- 1/2 system) if G can be embedded as a subgraph of the quantum hardware graph U. 
We refer this embedding problem as subgraph-embedding, to be defined formally in Section [3] In general, there 
are physical constraints on the hardware graph U. In particular, there is a degree-constraint in that each qubit 
can have at most a constant number of couplers dictated by hardware design. Therefore, besides the possible 
difficulty of the subgraph-embedding proble the graphs that can be solved on a given hardware graph U 
through subgraph-embedding must also be degree-bounded. Kaminsky et al. IITtI USl observed and proposed 
that one can embed G in [/ through ferromagnetic coupling "dummy vertices" to solve Maximum Independent 
Set (MIS) problem H of planar cubic graphs (regular graphs of degree-sjl on an adiabatic quantum computer. 
In particular, they proposed an n x n square lattice as a scalable hardware architecture on which all n / 3- vertex 
planar cubic graphs are embeddable. The notion of embedding here follows naturally from physicists' intuition 
that each logical qubit (corresponding to a vertex in the input graph) is mapped to a subtree of physical qubits 
(corresponding to vertices in the hardware graph) that are ferromagnetically coupled such that each subtree of 
physical qubits acts like a single logical qubit. For example, in Figure |2j the logical qubit 1 (in orange color) 
of the graph G is mapped to a subtree of physical qubits (labelled 1) of the square lattice. Informally, a minor- 
embedding Gemb of a graph G in the hardware graph [/ is a subgraph of U such that Gemb is an "expansion" 
of G by replacing each vertex of G with a (connected) subtree of U, or equivalently, G can be obtained from 
Gemb by contracting edges (same color in Figure |2l). In graph theory, G is called a (graph) minor of U (see for 
example |[T2l ). The minor-embedding will be formally defined in Section |3] (Remark: The embedding in ifTTlfTSl 
is a special case of minor-embedding, known as topological-minor embedding.) 

'Readers should be cautious not to confuse this embedding problem with the NP-compIete subgraph isomorphism problem, in which 
both graphs are unknown. However, in our case, the hardware graph is known. For example, if the hardware graph is a complete graph, 
then the embedding problem will be trivial. 

^MIS is a special case of QUBO and will be addressed in Section|5] 

^In their earlier paper 117J , it was said for graphs with degree at most 3, but the Ising Hamiltonian they used there was for regular 
graphs of degree-3. 




(3) 



iGV(G) 



iiGE(G) 
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Figure 2: Gemb(right) is a minor-embedding of G(left) in the square lattice U. Each vertex (called a logical qubit) 
of G is mapped to a (connected) subtree of (same color/label) vertices (called physical qubits) of U. G is called 
a (graph) minor of U. 

By reduction through minor-embedding, we mean that one can reduce the original Ising Hamiltonian on the 
input graph G to the embedded Ising Hamiltonian Ti,^"^^ on its minor-embedding Gemb> i-e-> the solution to the 
embedded Ising Hamiltonian gives rise to the solution to the original Ising Hamiltonian. The intuition suggests 
that the reduction will be correct provided that the ferromagnetic coupler strengths used are sufficiently strong 
(i.e., large negative number). However, how strong is "strong enough"? In lITTl [TSl . they do not address this 
question, i.e. what are the required strengths of these ferromagnetic couplers? In Section l4~n we will show that 
it is not difficult to give an upper bound for the ferromagnetic coupler strengths and thus explain the intuition. 
However, there are indications (T] that too strong ferromagnetic coupler strengths might slow down the adiabatic 
algorithm. Furthermore, an adiabatic quantum computer is an analog computer and analog parameters can only 
be set to a certain degree of precision (a condition much more stringent than the setting of digital parameters). 
Hence, the allowed values of coupler strengths are limited. Therefore, from the computational point of view, it 
is important to derive as small (in terms of magnitude) as possible sufficient condition for these ferromagnetic 
coupler strengths. Furthermore, what should the bias for physical qubits be? 

There are two components to the reduction: embedding and parameter setting. The embedding problem is 
to find a minor-embedding Gemb of a graph G in U. This problem is interdependent of the hardware graph 
design problem, which will be discussed in Section [6l The parameter setting problem is to set the corresponding 
parameters, qubit bias and coupler strengths, of the embedded Ising Hamiltonian. In this paper, we assume that 
the minor-embedding Gemb is given, and focus on the parameter setting problem (of the final Hamiltonian). Note 
that there are two aspects of efficiency of a reduction. One is how efficient one can reduce the original problem to 
the reduced problem. For example, here we are concerned how efficiently we can compute the minor-embedding 
and how efficiently we can compute the new parameters of the embedded Ising Hamiltonian. The other aspect 
concerns about the efficiency (in terms the running time) of the adiabatic algorithm for the reduced problem. In 
general, the latter depends on the former. For example, the running time of the adiabatic algorithm may depend 
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on a "good" embedding that is reduced to. According to the adiabatic theorem (see, e.g. ll22l ). the running 
time of the adiabatic algorithm depends on the minimum spectral gap (the difference between the two lowest 
energy levels) of the system Hamiltonian, which is defined by both initial Hamiltonian and final Hamiltonian. 
That is, the running time (and thus the efficiency of the reduction) will depend on both the initial Hamiltonian 
and the final Hamiltonian. In this paper, our focus is only on the final Hamiltonian, and therefore we are not 
able to address the running time of the adiabatic algorithm. Furthermore, the estimation of the minimum spectral 
gap of the (system) Hamiltonian is in general hard. Consequently, analytically analyzing the running time of an 
adiabatic algorithm is in general an open question. 

Finally, let us remark that there is another different approach based on perturbation theory by Oliveira & 
Terhal {2V\ for performing the reduction. In particular, they employed perturbative gadgets to reduce a 2-local 
(system) Hamiltonian to a 2-local Hamiltonian on a 2-D square lattice, and were able to show (as in the pioneering 
work 1191 ) that the minimum spectral gap (and thus the running time) of the system Hamiltonian is preserved (up 
to a polynomial factor) after the reduction, for any given initial Hamiltonian. However, besides the ineffective 
embedding, as pointed out in fTO\, the method is "unphysical" as it requires that each parameter grows with the 
system size. 

The rest of the paper is organized as follows. In Section |2l we recall the equivalences between the QUBO 
problem and the Ising problem. In Section [3l we introduce the minor-embedding definition and mention related 
work in graph theory. In Section |4j we derive the new parameters, namely values for the qubit bias and the 
sufficient condition for the ferromagnetic coupler strengths, for the embedded Ising Hamiltonian such that the 
original Ising problem can be correctly solved through the embedded Ising problem. In Section |5l we show the 
embedded Ising Hamiltonian for solving the (weighted) MIS problem. Finally, we conclude with several related 
algorithmic problems that need to be investigated in order to facilitate the design of efficient adiabatic algorithms 
and AQC architectures in Section |6] 



2 Equivalences Between QUBO and the Ising Problem 

In this section, we recall the equivalences between the problem of QUBO (maximization of 3^ in Eq. ([3])) and the 
Ising problem (minimization of £ in Eq. Notice that Xi = that is, Xj = 1 corresponds to Si = 1, and 
Xi = corresponds to Sj = —1. Using a change of variables, we have 

y{xi, . . . ,Xn) = - Ci — - Jij — - I I Jij — 2Ci j Si + JijSi.Sj j 

where nbr(i) =^ {j : ij £ E{G)}, the neighborhood of vertex i, for i G V(G). 

Therefore, Max3^(xi, . . . , x„) inEq. Q is equivalent to Min<S(si, . . . , s„) inEq. ^ where hi = X]jgnbr(i) 
2ci, or Ci = l/2(X]jgnbr(i) '^'■j ~ ^'^^ ^ ^ ^i^)- See Figure|4]for the correspondences between parameters in 
QUBO and the Ising model. 



3 Minor-Embedding 

Definition 3.1. Let U be a fixed hardware graph. Given G, the minor-embedding of G is defined by 

(p:G — >U 

such that 
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QUBO 



Ising Model 




(LJij -2q) 




-Ji 



1/2(1 J. - hi) 






Figure 3: The correspondences between the parameters in QUBO and the Ising Model. 



• each vertex in V(G) is mapped to a connected subtree Ti ofU; 

• there exists a map r : V(G) x V(G) — > V(C/) such that for each ij G E(G), there are corresponding 

e V(ri) and G \/{Tj) with ir{i,j)jr(j,i) ^ E(C/). 

Given G, if (j) exists, we say that G is embeddable in U. In graph theory, G is called a minor of U. When (p is 
clear from the context, we denote the minor-embedding (j){G) ofG by Gemb- 

Equivalently, one can think of a minor G of [/ as a graph that can be obtained from a subgraph of U by 
contracting edges. See Figure|2]for an example. 

In particular, there are two special cases of minor-embedding: 

• Subgraph-embedding: Each Tj consists of a single vertex in U. That is, G is isomorphic to G^mh (a 
subgraph of U). 

• Topological-minor-embedding: Each Tj is a chain (or path) of vertices in U. 

Minors are well-studied in graph theory, see for example |[T2ll . Given a. fixed graph G, there are algorithms 
that find a minor-embedding of G in [/ in polynomial time of size of U, from the pioneering 0(|V([/)p) time 
algorithm by Robertson and Seymour [23] to recent nearly linear time algorithm of B. Reed (not yet published). 
However, it is worthwhile to reiterate that these algorithms are for fixed G, and their running times are exponential 
in the size of G. Here the minor-embedding problem is to find a minor-embedding of G (for any given G) while 
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fixing U. To the best of our knowledge, the only known work related to our minor-embedding problem was by 
Kleinberg and Rubinfeld flOl, in which they showed that there is a randomized polynomial algorithm, based on 
a random walk, to find a minor-embedding in a given degree-bounded expander. 

Our embedding problem might appear similar to the embedding problem from parallel architecture studies. 
However, besides the different physical constraints for the design of architectures, the requirements are very 
different. In particular, in our embedding problem, we do not allow for load > 1, which is the maximum number 
of logical qubits mapped to a single physical qubit. Also, we require dilation, which is the maximum number 
of stretched edges (through other qubits), to be exactly 1. However, all of the existing research on embedding 
problems for parallel processors lfT3l . at least one of the conditions is violated (namely either load > 1 or dilation 
> 1). In this paper, our focus is on parameter setting such that the reduction is correct, and not on the minor- 
embedding algorithm and/or the related (minor-universal) hardware graph U design. These problems will be 
addressed in a subsequent paper. 



4 Parameter Requirement for the Embedded Ising Hamiltonian 

Let the QUBO problem, specified by y and G, be given as in Eq.©, and the corresponding £ in Eq.©. Suppose 
Gemb is a minor-embedding of G, and let £^'^^ be the embedded energy function associated with Gemb- 

£''"'^{si,...,Sn) = ^ Ksi+ ^ JijSiSj 

jeV{Gemb) «iGE{Ge,T,b) 

where \y{G,^k)\=N. 

Our goal is to find out the requirement for new parameters, h's and J's, such that one can solve the original 
Ising problem on G by solving the embedded Ising problem on its embedding Gemb> or equivalently, such that 
there is an one-one correspondence between the minimum of £ (and thus the maximum of y) and the minimum 
of £^"'''. 

As we mentioned in the introduction, the idea that we can solve the original Ising problem (i.e. finding the 
groundstate of the Ising Hamiltonian) on the input graph G by solving the new Ising problem on the embedded 
graph Gemb in U is that one can use ferromagnetic couplers to connect the physical qubits in each Tj of Gemb 
such that the subtree Tj will act as one logical qubit i of G. 



Notation. First, we recall that V(Gemb) = UieV(G) ^(T'O and E(Gemb) = UiGV(G) ^i'^i)^[JijeE(G) iT{i,j)3r{j,i)- 
We distinguish the edges within TjS from the edges corresponding to the edges in the original graph. Denote the 
latter by OE(Gemb), that is, OE(Gemb) = UijGE(G) V(ij)JrO,i)- ^^^^ ^^^^^ ^^8^^ '^^mb of Figure El corre- 
spond to OE(Gemb)-) For convenience, Onbr(ifc) =^ nbr(ifc) n OE(Gemb) and Fnbr(i,fc) =^ nbr(ifc) \ Onh<c{ik), 
ik G V(ri). Note j^(j-i) G 0nbr(4) ^ ij G E(G). 

The above intuition suggests that we use the same coupler strength for each original edge, i.e., J- ■ = 

l-T(i,j)jT(j,i) 

Jij for ir{ij)jT(j,i) ^ OE(Gemb)> and use ferromagnetic coupler strength, Ff{< 0), for each edge e G E(rj), and 
redistribute the bias hi of a logical qubit i to its physical qubits in Tj. That is, we choose h'-^ for physical qubit 
ik G y{Ti) such that J2i,&v{T,) = hi- 
Therefore, we have 

ieV(G) \ifcGV(r,) ipiq&E{T,) J iieE(G) 
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4.1 An Easy Upper Bound for the Ferromagnetic Coupler Strengths 



In this section, we derive an easy upper bound for the ferromagnetic coupler strengths. The derivation is based 
on the penalty or multiplier method, in which a constrained optimization is reduced to an unconstrained one by 
replacing the corresponding expression into the objective function as a penalty term with a large muhipUer0 For 
example, this was used in the reduction from general (higher-order) unconstrained binary optimization to the 
quadratic one (see ifSl). 

Notice that by construction (i.e., G is a minor of Ggmb and /i-^s are chosen such that ^jj.gv(Ti) Kk ~ 
minimization of <5 in Q is equivalent to 

min£^constrained(si, • • • ! STv) = E E ^ik^^k ~^ E '^*i'^*T(i,j) ■^JrO.i) 

i6V(G)ifeeV(T,) iieE(G) 

subject to 

SipSig = 1 for all ipig G E{Ti), i G V(G) 

where the condition Si^Si^ = 1 for all ipig S E(rj) is equivalent to requiring the spins of physical qubits that 
correspond to the same logical qubit i to be of the same sign. 
The corresponding unconstrained minimization is thus 

jeV(G) jfcev(7;;) \ ipi,eE(7-.) / ijeE(G) 

which is equivalent to solving Eq.© as X]ieV(G) X]jpjq6E(Ti) ^i"^ ^ constant. We are interested in how large 
(in terms of magnitude) F[^s are sufficient to guarantee that the solution of Eq.® gives the solution to Eq.©, 
and consequently to Eq.©. The result is stated in the following theorem. 

Theorem 4.1. Let fconstrained <^nd £' given as in Eq.^ and Eq.^. Suppose that for i G V(G), 

F[''<-l\hi\+ Yl foripig £E{Ti) (7) 

\ jenbr(i) / 

then we have s*^sl^ = Ifor all ipig S E(Tj), i G V(G), where {s^^, . . . , s^) = arg min£^'. Consequently, 

mill ^constrained (Sl, ■■■,Sn) = min<S'(si, ...,Sn)- 



Proof. Let (sj^ . . . , s*^) = arg min£^'. Suppose on the contrary. That is, there exists fpfg G E{Tf ), for some 
1 < f < n such that s*f s*f = — 1. Then we have 

— •' — Jp Iq 

\/fceV(T^) jGnbr(/) / 

"'This is closely related to the discrete Lagragian Multiplier method. However, the Lagragian multipliers are treated as unknown or 
iteratively solved with the Lagragian function. 
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where 

ieV(G),i^/ifceV(T,) \ ipigeE(T,) / ijeE{G),iytf 

Thus, if Fs satisfy Eq.(??), we have 

/fceV(T^) jGnbr{/) 

which is at least min£^'(si, . . . , sn) = £'{s1, . . . , s*j^), a contradiction. Consequently, we can conclude that 
s* s* =1 for all ipiq € E(Tj), i € V(G), where (s^^, . . . , s^) = arg uiinS', implying the claimed equality of 
the minima. □ 



4.2 A Tighter Bound for the Ferromagnetic Coupler Strengths 

As we mentioned Section \T\ it is desired to obtain tighter bounds for the ferromagnetic coupler strengths. In 
this section, we show that with a more careful analysis, we can reduce the bound by setting the qubit bias (h') 
appropriately. 

Fori G V(G),let 

Ci — ^ ^ I'^ijl I ^ j I • 
jenbr(i) 

Observe that if Ci < 0, that is, X]jgnbr(j) ^ 1^*1- Then we have s* = —1 for hi > 0, and s* = +1 for 
hi < 0, where (s^, . . . , s*) = arg mm£. Therefore, WLOG, for the rest of the paper, we will assume that 
Ci>0. 

Theorem 4.2. Let £ (y resp.) and G be given as in Eq.^ (Eq.^ resp.}. Suppose Ggmb is a minor-embedding 
of G, and let £^'^^ be the energy of the embedded Ising Hamihonian given in Eq.@. Then for all i € V((S'), if 



and 



■ ,u\\ ^jV(, rtSOnbrfe) \Ji \ - Gi/k Ik IS one of the k leaves ofTi 
/i- = signihi) { ■'^''J'^i I T I 1 

«fc ^ \ \ y ^n„hr(i \\Jii\ otherwise 



Ff < --\^Gi for all E{T,), 



where Cj > (defined above) and sign{hi) = +1( — 1 resp.) if hi > (hi < resp.), we have s*^s*^ = 
1 for all ipiq € E(Tj), i G V(G), where {s*_^, . . . , s^) = arg minf^"^'^. Consequently, there is one-one corre- 
spondence between arg min£'^'^'^ and arg minf (and thus arg maxj^j. 

Furthermore, if we set Ff = — ^-^f-^Ci — gi/2, for some Qi > 0, then the spectral gap (which is the difference 
between the two lowest energy levels) of the embedded Ising Hamihonian will be the minimum o/minjgv(G) 9i 
and the spectral gap of the original Ising Hamiltonian. 

See Figure |4] for an example of the corresponding parameters in the embedded Ising Hamiltonian (for the 
topological-minor-embedding case). 

In the following, we first explain the main idea behind, followed by the formal proof. For the illustration 
purpose, suppose hi > 0, consider the simplest case in which all but one leaf is —1. Now, consider the energy 
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change if the leaf is flipping from +1 to —1. Our goal is to set Fs as small (in terms of the magnitude) as 
possible, such that the energy change is at least greater than zero. Note that the energy change A£i^ = 2(/i^^ + 
Ejr(i,i)GOnbrfe) ^Mj.) Jij " ^) ^hcrc 6j^^^^^^ depends on the sign of spin s*^^^^^. We would like to bound A£i^ 
without knowing the signs of Sj s. Observe that if Js are all positive, then the worst case we have A£^ 



> 



2(M,-E,,(„,)eOnbrfe) J^J-F). And in general A^,, > 2(/i^^-E,^^^ ^.^^o^^rfe) 4+Sj.(,„)e0nbrfe) J^J-F), 
where J,rt^ > and J-^ < 0. Consequently, setting F <h^^- Y^j^^^ j)GOnbr{ifc) \ \ would imply A£i^ > 0. One 
can then extend this argument to a segment that needs to be flipped, as illustrated in Figure [5] 

Proof. It is easy to check that J2ii,G\/{Ti) Kk ~ sign{hi)\hi\ = hi. 

We wiU first prove the theorem when Gemb is a topological-minor-embedding of G. That is, each Tj is a 
chain (ii, • • • > ^tj connected by consecutive vertices. Thus, the corresponding embedded energy function is 
given by 

/ U ti-l \ 

jGV{G) \A:=1 k=l / jieE(G) 

In this case, we have li = 2. We prove by contradiction. Suppose NOT. That is, there exists i e V(G), and 
1 < k < ti such that s*^s*^^_^ = — 1. We distinguish two cases based on the value of hi. 

For hi > 0, let p be the smallest index such that s* = — 1 and s* = +1, and let q be the smallest index 
such that s* = +1 and s* = -1. That is, we have s* = . . . = s* = -1, s* = . . . = s* = +1. Note 

tq 1 ''1 ''p—1 *p f-q 

we have at least either p > 1 or g < n. (Otherwise they are all +1). Then we claim that by flipping the segment 
{ip, . . . ,ig) from +ls to —Is, the energy decreases, conti"adicting to the optimality. (Notice that the Fs within 
the segment do not change.) The energy change equals to 

> 2(^2^,-^2 E l-^^i I - ^i"''^" - ^f'^^^ ] ^^^t*^ convention that F^^ = F^'^^^"^ 

\k=p fc=pjV(j ,j)GOnbr(jfc) / 

> 2(-Ci/2 - F^''^'^^'') (where k=pifp>l,k = q otherwise) 

The argument for hi < Ois similar except that we will flip from —1 to +1 instead, and 

^^emb ^ ^emb(_ . . , +1, . . . , +1, _1, . . . , _1, +1, . . .) _ £er.b^_ . . , +1, . . . , +1, +1, . . . , +1, +1, . . .) 

^ V ' 

> 2(-Ci/2 - fI'^''^^'^) (where k = p if p>l,k = q otherwise) 

Therefore, in both cases, if Ff < -Ci/2, for all 1 < p < U, we have A£^"''° > 0, contradicting to the 
optimality. Hence all s*^ must be of the same sign. 

For the general minor-embedding when TjS are trees, we can similarly argue that for hi > 0, if all but one 
leaf is positive (+1), then flipping them (from +ls to —Is) will decease the energy provided that Fj < — ^V^Cj. 
Similarly, one can argue for the case when hi < 0. 
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Notice that if we set Fi = —^-j^Ci — gi/2, then the energy change lS.£^"^^ > 5j- In this case, the spectral 
gap of the embedded Ising Hamiltonian will be the minimum of minjgv(G) 9i ^^i^ spectral gap of the original 
Ising Hamiltonian. □ 

Remark. One can generally set h'^^ = sign{hi){^- ^ .jgOnbr(ifc) ~ ^ik) where Cij.s are chosen such that 
X^ifceV(Tj) ~ C'j(= Sjgnbr(j) \ ~ l^d)' ^^^^ ^^c Fs accordingly. This flexibility is useful when the 
precision for parameters is limited. (In the above theorem, we give an upper bound for the ferromagnetic coupler 
strengths. A natural question is how tight our bound is. Can one get an better bound without solving the original 
problem?) 

5 Weighted Maximum Independent Set (WMIS) Problem 

First, we formulate WMIS problem as a special case of QUBO in Section [STTl Then, we apply Theorem 14.21 to 
set parameters for the embedded Ising Hamiltonian for the corresponding MIS problem in Section l5!2l 

5.1 Formulate WMIS Problem As a Special Case of QUBO 

Weighted MIS (WMIS). Given an undirected vertex-weighted graph G. Let V(G) = {1,2,..., n}, let a G 
M"*" be the weight of vertex i. WMIS seeks to find a wmis(G) = S QV such that S is independent and the total 
weight of S (= X]ie5 ^j) is maximized. 

Tlieorem 5.1. If Jij > minjcj, Cj}/or all ij G E(G), then the maximum value of 

y(x\, . . . , Xn) — ^ ^ CiXi ^ ^ JijXiXj 
iGV(G) «jeE(G) 

is the total weighted of the WMIS. In particular if Jij > min{cj, Cj} for all ij G E(G), then wmis(G) = {i G 
V(G) : X* = 1}, where {x\, . . .,x^) = arg max(3.^^„,^^^)g{o,i}" Vi^i, ■ ■ ■,Xn)- 

Proof of Tlieorem 15.11 Let {xl, . . . = arg max^^.^^ ,j,^-)g|Q ]^|„ y{xi, . . . ,x„). Denote S* = {i e V(G) : 
X* = 1}. We'll prove that if Jij > minjcj, Cj} for all ij G E(G), then S* is an independent set. 

Suppose on the contrary, that is, there exists an edge yz in the subgraph induced by S*. WLOG, as- 
sume that Cy < Cz- Consider removing y from S*. Let S' = S* \ {y}. The weight change equals to 
-Cy + J2jenbr{y)ns* -^yJ - ~^ + "^yz > 0' contradicting to the optimahty of S* . □ 

The above theorem is a generalization of the known fact for unweighted case of MIS (see fS] and references 
therein). For the unweighted case of MIS, Cj = 1 for all i G V(G). Thus, it is sufficient to choose Jij = 1 + e 
for all ij G E(G) for some e > 0. Accordingly, the corresponding energy function of the Ising Model for MIS is 
£{si,...,Sn) = E»gv(G) (degi(l + e) - 2) s, - EijGE(G)(l + 

Remark. If we choose J to be exactly 1 instead (as in lUTl [T8l|), then we can only guarantee the size of the 
maximum independent set, but the returned set is not necessarily independent. For instance, when G = i^^4, any 
(adjacent) two vertices also has the minimum energy of one. 

Note that we thus can conclude that the Ising problem is NP-hard because WMIS is NP-hard. Indeed, Bara- 
hona [71 showed the NP-hardness of a special Ising problem through the reduction of MIS problem on cubic 
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graph (which remains NP-hard). Notice that for WMIS on a planar graph, there is a PTAS algorithm Q. Re- 
cently, for the Ising problem on a planar graph, Bansal et al. [6| applied the same technique to obtain a PTAS 
algorithm in 0(n2^^/'^) where (1 — 5) is the approximation ratio. As another side note, conversely (to the fact 
that WMIS is a special case of QUBO), Boros et al (see e.g. JH) has shown that a problem of QUBO on G can 
also be converted to a WMIS but on a different graph. 

5.2 Embedded Ising Hamiltonian for Solving MIS 

Let G be the graph of MIS problem and Gemb be the topological-minor-embedding of G in U. For the unweighted 
MIS, we have Ci = X]jGnbr(i) — 2- Therefore, according to Theorem l4.2[ it suffices to set Fi < —Gi/2 = 

—1. That is, the embedded Ising Hamiltonian for solving unweighted MIS is given as in Eq. (H), with 

1/ J dii^J — 1 k = l,ti 

~\ d^,J l<k<ti 

where di^. = |Onbr(zfc)|, and Jij = J > 1, F^'^^^^'^ = F < —1. In particular, for degree-3 hardware graph U, 
by setting J = 1 + e, and F = — (1 + e), for some e > 0, there are only 6 different parameter values, namely, 
{-(1 + e), 0, e, 1 + e, 1 + 2e, 1 + 3e} in Eq. dUl. 

6 Discussion 

In this paper, we introduce minor-embedding in AQC. In particular, we show that the NP-hard QUBO problem 
can be solved using an adiabatic quantum computer that implements Ising spin- 1/2 Hamiltonians, through minor- 
embedding reduction. There are two components to this reduction: embedding and parameter setting. Given a 
minor-embedding, we show how to derive the values for the corresponding parameters, in particular, a good 
upper bound (in terms of magnitude) for the ferromagnetic coupler strengths, of the embedded Ising Hamiltonian 
such that there is one-one correspondence between the groundstate of the original Ising Hamiltonian and the one 
of the embedded Ising Hamiltonian. 

There are many algorithmic problems related to minor-embedding in AQC that remain to be addressed. In 
particular, the problems relate to the efficiency of the reduction. These problems in turn relate to the running 
time or complexity of quantum adiabatic algorithms. Recall that according to the adiabatic theorem, the running 
time of an adiabatic algorithm depends on the minimum spectral gap of the system Hamiltonian, which however 
might be as hard as solving the original problem. Despite several serious investigations, the power of AQC 
remains an open question ll24ll25ll22l[T6l [Tl. How does the embedding reduction effect the time complexity of an 
adiabatic algorithm? In order to address this question, one will also need to specify the initial Hamiltonian. In [2l|, 
we show that for some special cases, how the embeddings, parameters, and initial Hamiltonians can effect the 
minimum spectral gaps. The effect of the embedding and its consequential initial Hamiltonian on the complexity 
of adiabatic algorithms remains to be investigated. In the following we state several main problems that need to 
investigate in order to facilitate the design of adiabatic algorithms and AQC architectures. Partial results to these 
problems will appear in our subsequent papers. 

PI. Measurement for the minor-embedding. Define a measure for the minor-embedding such that a good 
minor-embedding corresponds to a reduced problem that admits an efficient adiabatic algorithm. 
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P2. Embedding-dependent initial Hamiltonian. Design an embedding-dependent initial Hamiltonian for a 
given minor-embedding such that the adiabatic algorithm for the reduced problem is at least as efficient as the 
adiabatic algorithm (with the best possible initial Hamiltonian) for the original problem. 

P3. Hardware graph design. Given a family T of graphs (which consists of classically hard instances), the 
problem is to design a hardware graph (called a T -minor-universal graph) which is as small as possible (in terms 
of total number of vertices and edges) such that 

• all known physical constraints are satisfied; 

• all graphs in T are embeddable; 

• a good embedding of each graph in T can be efficiently computed. 
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Appendix 



A Energy Function of the Ising Hamiltonian 

In this section, we show that the energy function of the Ising Hamiltonian in Eq. ^ is given in Eq. That 
is, replace each uf by the variable Sj G {— +1}- The equivalence follows from two basics: (1) Eigenvalues 

and eigenvectors of a^: a^\0) = (+1)|0), a'\l) = (-l)Il), where |0) =^ [l,0]t, |1) =^ [0,1]^. (2) Tensor 
product property: (^4 (g) B){\u) (g) \v)) = A\u) ® B\v). More precisely, for \z) = \zi) . . . (g) Zi G {0, 1}, 
i = 1, . . . ,n. 

^finalk) = Yl ^i^t+ Y Jij(^t<^'j\\zi) ® ■■■®\Zn) 

\iGV(G) iieE(G) / 

= Y h\zi)® ...®{-iY^\z,) ...®\zn)+ Y JijW) ® ■■■'^i-^y"\^i) ■■■®i-'^y'\^i) 

^iGV{G) jj6E(G) J 

Therefore, the energy function of Hf\na\ is 

£{zi,...,zn)= Y h^(-^y'+ E M-^n-^p 

iGV{G) iieE(G) 

where Zi G {0, 1}, for i = 1, . . . , n. Replace (—1)^' by Si € {+1, —1}, we thus have Eq.©. 
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(generalJ) (J>0) 

Figure 4: The corresponding parameters in the embedded Ising Model. Bottom right, for J > 0; bottom left, for 
general J: Jn > 0, Ji2 < 0, J^a < 0, and = Q - Ejenbr(*),j,,<o ^ij- 
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{c)h > {d)h < 



Figure 5: (a) For h > and positive Js (denoted by J"*"), the energy change AS > 2{h — — — F) when 
flipping the leaf's spin from +1 to —1. (6) For /i > 0, /S.8 > 2{h — + J~ — F) when flipping the leaf's spin 
from +1 to —1, where J+ > 0, < 0. (c) For /i > 0, flip a segment from +1 to —1. (d) For /i < 0, flip a 
segment from —1 to +1. 
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